Multi-gap discrete vector solitons 
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We analyze nonlinear collective effects in periodic systems with multi-gap transmission spectra 
such as light in waveguide arrays or Bose-Einstein condensates in optical lattices. We demonstrate 
that the inter-band interactions in nonlinear periodic gratings can be efficiently managed by con- 
trolling their geometry, and predict novel types of discrete vector solitons supported by nonlinear 
coupling between different bandgaps and study their stability. 
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Periodic structures are common in nature, with the 
crystalline lattice being the most familiar example. One 
of the important common features of such systems is the 
existence of frequency gaps in the transmission spectra 
which can dramatically affect both propagation and lo- 
calization of waves. Moreover, the modern technology 
allows creating different structures with an artificial pe- 
riodicity, and the recent examples are photonic crystals, 
which can control propagation and emission of electro- 
magnetic waves fl| , and optical lattices, which are used to 
trap and manipulate atomic Bose-Einstein condensates 
(BECs) 0. An unprecedented level of control over such 
engineered structures can be realized by tailoring both 
location and width of multiple band-gaps with additional 
modulation of the structure parameters. For example, it 
has been shown that atomic BEG can demonstrate a rich 
variety of phase transitions in optical superlattices 
and the reduced-symmetry photonic crystals allow self- 
localization of waves in mini-gaps . 

The response of many systems becomes nonlinear at 
higher energy densities. This phenomenon may have var- 
ious physical origins, such as the charge recombination 
in biased photorefractive crystals 5j , excitation of higher 
energy levels in semiconductors 6] , or atomic interaction 
in BEG In periodic media, nonlincarity produces a 
shift of the band-gap spectrum, and this physical mecha- 
nism is responsible for a number of remarkable effects, in- 
cluding the formation of gap solitons 8] . Such nonlinear 
localized modes can be excited within multiple spectral 
gaps of a periodic structure, as was first demonstrated 
experimentally for temporal optical pulses in fiber Bragg 
gratings Since the pulses extend over hundreds of 
grating periods, the gap-soliton dynamics is usually de- 
scribed by avera ged coupled-mode equations with con- 
stant coefficients [lO|. In contrast, spatial optical beams 
in waveguide arrays and matter waves in BEG can span 
over a few periods of the structure. Under such con- 
ditions the wave profiles are essentially discretized by 
the underlying periodic structure, strongly affecting the 
properties of discrete solitons [lllll^. 

Spatial optical solitons associated with the first spec- 
tral band of the multi-band transmission spectrum have 
been extensively studied both theoretically and exper- 
imentally in arrays of coupled optical waveguides [llj. 



Very recently, spatial gap solitons localized in higher- 
order bands have been observed as well J^]. Similar 
observation of matter-wave gap solitons in BEG loaded 
into an optical lattice is expected soon as well. How- 
ever, the interaction properties of localized modes which 
belong to different gaps are not known. The existence 
of gap solitons is directly related to the band-gap spec- 
trum, and the latter can be fine-tuned in superlattices. 
In this Letter, we reveal, for the first time to our knowl- 
edge, the fundamental links between periodic modulation 
of the medium parameters and nonlinear wave coupling 
between different gaps, and also predict the existence of 
multi-gap discrete vector solitons with nontrivial symme- 
try and stability properties. We believe that our results 
can stimulate and guide the future experiments on dis- 
crete gap solitons in optics and matter-wave physics. 

Self-action and interaction of optical beams in a one- 
dimensional periodic structure of coupled optical waveg- 
uides with the normalized refractive index profile V{x) — 
V[x + h) can be described by a set of coupled nonhnear 
Schrodinger equations. 



dz dx^ 



V{x)i}} + G(V')V' = 0, 



(1) 



where, in the case of optical waveguides, x and z are the 
transverse and longitudinal coordinates, respectively, h is 
the spatial period, ij) = {ip'^^\ip'^'^\ ■•■) ^re the normalized 
electric field envelopes of several co-propagating beams 
having different polarizations or detuned frequencies. It 
is assumed that the beams interact incoherently through 
the Kerr-type nonlinear change of the refractive index, 
Gm = Y.j rmj#'^' where Tm=j and Tm^j are the self- 
and cross-phase modulation coefficients, respectively. 

We note that the model is equivalent to a system 
of the coupled Gross-Pitaevskii equations describing the 
dynamics of multi-component BEG in ID optical lattice, 
where -0^'") is the mean- field wave function for atoms in 
the m-th quantum state, z stands for time, V{x) is the 
periodic potential of an optical lattice, and G is the effec- 
tive mean-field nonlincarity which appears due to the s- 
wave atom interaction. Although below we use the termi- 
nology from guided wave optics where the rapid progress 
in the experimental study of spatial optical solitons is ob- 
served, our results are equally applicable to the nonlinear 
dynamics of BEG in optical lattices. 



Linear wave propagation through a periodic structure 
can be entirely described by the Floquet-Bloch spec- 
trum of the eigenmode solutions of Eqs. (Q) in the form 
"0 = iph 6x^(1 f3z + iKfjx/h), where K}, is the normalized 
Bloch-wave number and (3 is the propagation constant. 
The Bloch wave amplitudes decay exponentially when 
lmKt,{/3) ^ 0, and this condition defines the location of 
gaps in the transmission spectrum. At the gap edges, 

Nonlinearity manifests itself through an effective 
change of the optical refractive index which results in 
a local shift of the bands and gaps. This physical mech- 
anism is responsible for the formation of gap solitons. 
When the band shifts are small, we can seek solutions of 
Eqs. © near the gap edges {(3 = Pm) in the form of mod- 
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<^('")4"' exp(i/3„z + 
iK^"^' x/h), and obtain a system of coupled nonlinear 
Schrodinger equations for the slowly varying envelopes. 



ulated Bloch waves [ij, V^"^ 
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^.if^WlV"^ =0.(2) 



Here 



are the effective diffrac- 



tion coefficients, and 7„jj = /J* I'^I^VI"'^ \'^dx are the non- 
linear coupling coefficients, where we assume the normal- 
ization /g'' IV'I™'' P'^^; = 1. It can be demonstrated 
that the diffraction coefficients D^™) are positive near 
the lower gap edges, and negative at the upper edges 
and, therefore, both bright and dark solitons can co-exist 
in the nonlinear media with either self-focusing or self- 
defocusing nonlinearities. Moreover, it immediately fol- 
lows from Eqs. ^ that each soliton can support multiple 
guided modes in other band gaps, all such modes can be 
coupled together to form multi-gap vector solitons. 

The simplified model (0) predicts stability of bright 
vector solitons 0| when all r„y have the same sign. 
However, the applicability of Eqs. Q is limited to small- 
amplitude solitons in the vicinity of the gap edges. Even 
in this regime, the important effect of the inter-band res- 
onances, which can lead to the soliton instabilities |l7l| . 
is not taken into account by the envelope approxima- 
tion (|2Jl. As the input power grows, the soliton width 
decreases and becomes comparable to the spatial period 
of the structure. This suggests that nonlinear properties 
should depend on the actual profile of the periodic poten- 
tial V{x). Modulated periodic structures or superlattices 
can then become an important tool to engineer both lin- 
ear and nonlinear properties of the Bloch waves and gap 
solitons. 

As an example, we consider a binary superlattice where 
the effective periodic potential is composed of two types 
(A and B) of separated individual potential wells, V{x) — 
X]„[Vk(a: 4- nh) Vb{x -\- nh)]. Such superlattices can 
be produced by etching waveguides on top of a AlGaAs 
substrate Il8l |. or induced dynamically by two overlap- 
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FIG. 1: Examples of a binary superlattice which can be cre- 
ated by (a) an array of two types of coupled waveguides, or (b) 
an optical superlattice induced by two overlapping mutually 
incoherent interference patterns. 
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FIG. 2: (a) Characteristic dependence of the Bloch wave 
number (Ki,) on the propagation constant /?. Gray shad- 
ings mark the transmission bands, (b) Dependence of the 
normalized self- and cross-phase modulation nonlinear coeffi- 
cients (7™^ = laf^^a'^'p -\- |6('"'6(-''p) between the gap edges 
/3i = P+ and (32 — ~I3- vs. the parameter p. The values of 
/3± correspond to the plot (a) by a proper choice of k. In- 
sets show possible symmetries of superlattices corresponding 
to different parameter values, but the same linear dispersion. 



ping mutually incoherent interference patterns in a pho- 
torefractive medium [l9j, see Figs. Ha,b). In order to 
analyze the properties of nonlinear waves of such su- 
perlattices, we employ the tight-binding approximation. 
This approach allows us to describe correctly the first two 
spectral bands, and the Rowland ghost gap ^| which is 
defined by a difference between the A- and B-type lat- 
tice cites. We present the total field as a superposition 
of the guided modes supported by individual potential 
wells, %l){x,z) = I]„[A„(z)V'A(a;) -I- B„(z)V'B(a;)], where 
A„ and Bn are (yet unknown) mode amplitudes. Finally, 
we derive a system of coupled discrete equations for the 
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normalized amplitudes an and b„, 

t-iCl^ -I 11^ 

i— h pa„ + b„_i + Kb„ + Xa|a„| a„ = 0, 

i-^- pb„ + Ka„ + K a„+i + Xb|b„| b„ = 0. 

Here, the key characteristics of the binary superlattice 
are defined by free parameters: p is proportional to the 
the detuning between the propagation constants of the 
A and B-type guided modes, n characterizes the relative 
coupling strength between the neighboring wells on the 
right- and left-hand sides, and Xa,b are the normalized 
nonlinear coefficients. 

According to Eqs. Q, the linear Bloch-wave dispersion 
is defined as Kb = cos~^(— 77/2), where = + -|- 

— (3. The transmission bands correspond to real Kb, 
and they appear when /3_ < < /?+, where /3± = 
(k^ + k^^-|-P^±2)^/^. a characteristic dispersion relation 
and the corresponding band-gap structure are presented 
in Fig. a). The upper gap at /3 > /?+ is due to the effect 
of total internal reflection (IR), whereas additional gaps 
appear due do the resonant Bragg reflection (BR). 

It is well known that the material dispersion can be 
completely compensated by the geometrical dispersion 
in optical fibers [i^ . More recently, diffraction manage- 
ment was realized in a periodic waveguide array struc- 
ture |23 |. The natural open question is whether it is 
possible to control nonlinear coupling between the gaps 
by appropriate design of periodic structures. In order to 
answer this fundamental question, we study the depen- 
dence of the nonlinear coupling coefficients on the super- 
lattice parameter p, while preserving exactly the same 
linear dispersion of Bloch-waves by a proper choice of k. 
Our results are presented in Fig. |2Ib), and they uncover 
the remarkable feature: nonlinear inter-hand interaction 
coefficients strongly depend on the symmetry of the peri- 
odic structure, and this relation cannot be fully charac- 
terized just by the linear Bloch-wave dispersion. Thus, 
by changing the superlattice parameters it is possible to 
selectively enhance or suppress inter-band interaction. 

To be specific, we consider the superlattice with sym- 
metric inter-site coupling {n — 1) created in a self- 
focusing medium. Such a lattice can support two funda- 
mental types of one-component bright solitons centered 
at either A or B cites, and these solitons exist in both 
the IR and BR gaps described by the model (O . We find 
that the powers of A- and B-type solitons become signif- 
icantly different away from the band edges, see Fig. Of a). 
The solitons of type-B in the IR gap are always unstable 
with respect to a translational shift (symmetry breaking) , 
however the stability is reversed for discrete gap solitons 
in the first BR gap [left part of Fig. OJa)] where type- 
A solitons become unstable. Additionally, the discrete 
gap solitons become oscillatory unstable above a critical 
power due to (i) internal resonance within the gap, first 
discovered for the fiber Bragg solitons , and (ii) inter- 




FIG. 3: (a) Power vs. propagation constant for discrete soli- 
tons centered at A (black) and B (gray) lattice sites. Solid — 
stable, dashed — unstable, and dotted — oscillatory unsta- 
ble modes. (b,c) Eigenvalues of the guided modes supported 
by the discrete solitons localized in the complimentary gap. 
Insets show characteristic profiles of solitons and their guided 
modes defined as U2n = an and U2n+\ = bn. The array pa- 
rameters are p = 0.75, k = 1, and Xafi = 1- 

band resonances, first found for nonlinear defect modes 
in a layered medium jlTj. 

The mutual trapping of the modes localized in dif- 
ferent gaps and the physics of multi-gap vector soli- 
tons can be understood in terms of the soliton-induced 
waveguides ■ Therefore, the effect of discreteness on 
the inter-gap coupling can be captured by studying the 
guided modes supported by a scalar soliton in other gaps: 
the larger is the eigenvalue shift from the band edge, 
the stronger is the interaction. In Figs. |3Ib,c), we plot 
the eigenvalues of the guided modes supported by the 
BR (gap) and IR solitons and observe two remarkable ef- 
fects which cannot be captured by the simplified envelope 
approximation First, the strength of the inter-gap 
coupling depends strongly on the soliton symmetry. In- 
deed, the type B soliton always creates a stronger effec- 
tive waveguide, despite the fact that the soliton power 
in the BR gap is smaller compared to the type-A soli- 
tons. Second, the nonlinear inter-gap coupling decreases 
for strongly localized discrete soliton in the IR regime, as 
follows from the non-monotonic dependence of the gap- 
mode eigenvalues shown in Fig. Efc) . 

The eigenvalues of the linear guided modes define the 
point where a multi-gap vector soliton bifurcate from 
their scalar counterparts. Initially the amplitude of the 
guided mode is very small, but it increases away from the 
bifurcation point, and the mode interacts with the soliton 
waveguide creating a coupled inter-gap state, see Fig. ^ 
In the vicinity of the bifurcation point, the soliton sym- 
metry and stability are defined by the large-amplitude 
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FIG. 4: (a,b) Powers in the IR (a) and BR (b) components 
of the muhi-gap discrete vector solitons vs. propagation con- 
stant /3i with /32 = 0.5 — /3i/3 for the families with symmet- 
ric (AA) and asymmetric (AB) profiles. Dashed lines mark 
solutions exhibiting symmetry-breaking instability. Bottom: 
Characteristic profiles of the vector solitons composed of the 
components localized in two different gaps: IR - unstaggered, 
and BR - staggered; propagation constants correspond to the 
points marked as I, II, and III in the upper plots. The array 
parameters match Fig.|21 



soliton component. For example, the AA-type discrete 
vector soliton shown in Fig. 01X11) is stable because the 
powerful A-type mode in the IR gap suppresses insta- 
bility of the second component. However, as the power 
in the second component grows, the soliton properties 
change dramatically: (i) AA state becomes unstable, 
and at the same time (ii) a stable AB-type asymmet- 
ric vector soliton emerges, see the modes II and I in 
Figs. Elbottom), respectively. These complex existence 
and stability properties underline a nontrivial nature of 
nonlinear inter-gap coupling between the localized com- 
ponents with different symmetries. 

We notice that such a simple way to engineer nonlinear 
coupling in the superlattices can lead to novel effects in 
the soliton collisions. Both coherent interaction of soli- 
tons from one band and vector solitons with the com- 
ponents from different bands can be controlled by engi- 
neering the superlattice parameters, thus leading to novel 
features in the soliton switching and steering. These in- 
teresting effects will be discussed elsewhere. 

In conclusion, we have studied nonlinear coupling and 



wave localization in periodic systems with multi-gap 
transmission spectra. We have predicted the existence 
of novel types of multi-gap vector solitons and studied 
their stability. Using the example of a binary waveguide 
array, we have demonstrated the basic concepts of the 
engineering of nonlinear inter-band interaction in such 
structures, which in turn determine the key soliton prop- 
erties. 

The authors acknowledge a support of the Australian 
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trovskaya. 
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